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ABSTRACT 

Context. The Optical Monitoring Camera (OMC) onboard INTEGRAL provides photometry in the Johnson V-band. With an aperture 
of 50 mm and a field of view of 5°x5°, OMC is able to detect optical sources brighter than V ~ 18, from a previously selected list of 
potential targets of interest. After more than nine years of observations, the OMC database contains light curves for more than 70 000 
sources (with more than 50 photometric points each). 

Aims. The objectives of this work have been to characterize the potential variability of the objects monitored by OMC, to identify 

periodic sources and to compute their periods, taking advantage of the stability and long monitoring time of the OMC. 

Methods. To detect potential variability, we have performed a chi-squared test, finding 5263 variable sources out of an initial sample 

of 6071 objects with good photometric quality and more than 300 data points each. We have studied the potential periodicity of these 

sources using a method based on the phase dispersion minimization technique, optimized to handle light curves with very dilferent 

shapes. 

Results. In this first catalogue of variable sources observed by OMC, we provide for each object the median of the visual magnitude, 
the magnitude at maximum and minimum brightness in the light curve during the window of observations, the period, when found, as 
well as the complete intrinsic and period-folded light curves, together with some additional ancillary data. 
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1. Introduction 

The study of variable stars is one of the principal areas of astro- 
nomical research. The beginning of their study dates back to the 
16th century when the first variable star, a supernova, was iden- 
tified, and shoitly thereafter telescopes started to be used for as- 
tronomy. Since then, an increasing number of variable stars were 
discovered and observed but it was with the upcoming of first 
photography and photoelectric photometry and later of the CCD 
technology and multi-wavelength astronomy that the number of 
known variable stars increased very rapidly. In 1948 the first 
edition of the General Catalogue of Variable Stars (GCVS) was 
published containing 10 820 stars. The GCVS has been updated 
regularly until its last release ( Samus et al.|2012) . The number 
of known variable stars has increased up to around 200 000 ob- 
jects compiled in th e Variable Star Index (VSX; |Watson||2006[ 
[Watson et al. 2012[ l. The VSX is maintained by the American 
Association of Variable Stars Observers (AAVSO), and contains 
all objects from the GCVS and the catalogues of the Noithern 
Sky VariabiHty Survey (NSVS; [Wozniak et al.|2004| l, as well as 
data resulting from other automated surveys, like the All Sky 
Automated Survey (ASAS; Pojmanski 1997), or the catalogues 
compiled from the Optical Gravitational Lensing Experiment 
(OGLE; |Udalski et al.|1992i ) observations; the Miras and eclips- 
ing binaries from OGLE-II data or the catalogue of variable stars 
from OGLE-III. 



The advent of space astronomy has provided additional high 
quality data for thousands of objects. The most significant im- 
pact has been produced by the missions specifically dedicated 
to obtain photometric data of high accuracy, in some cases dur- 
ing very long periods of time. While optimized for astrometric 
studies, the Hipparcos/Tycho mission also produced high accu- 
racy photometric time series ( Ferryman et al.|1 997, Hoeg et al. 
1997). Even higher photometric accuracy is provided by ded- 
icated missions, aimed to study asteroseismology or to search 
for exoplanetary transits, like the Canada's Microvariability and 
Oscillations of STars (MOST; 'Walke r et aLl|2003]) the W ide- 
field Infrared Explorer (WIRE; Br untt & Buzasi] |2006| , the 
convection, Rotation and planetary Transits (CoRoT) mission 
( [Auvergne efal. 2009,) or Kepler ( |Koch et al!][2010 j, which is 
yielding the highest accuracy Ught curves available up to now. 

Furthermore, spacecraft systems or secondary instruments 
on the satellites can provide very useful photometric informa- 
tion as well. Zwintz et al. ( 2000 1 performed a detailed analysis 
of more than 4500 data sets from the Hubble Space Telescope 
(HST) guide star data, finding 20 stars to be variable. In a similar 
way, the observations performed by the spacecraft tracking cam- 
eras of Chandra have allowed to produce the Chandra Variable 
Guide Star Catalogue (VGUIDE) containing data of 827 stars 
for which accurate, long term photometry was available (Nichols 
leFaL 2010). 



* Available in electronic form at the CDS via anonymous ftp to 
cdsarc.u-strasbg.fr (130.79.128.5) or via http://cdsweb.u-strasbg.fr/cgi- 
bin/qcat?J/A+A/ and at at http://sdc.cab.inta-csic.es/omc/\>ar/ 



The International Gamma-Ray Astrophysics Laboratory 
(INTEGRAL ) mission, launched October 17th, 2002, ( |WinkIer| 
et al. 2003 ) includes a small optical telescope, the Optical 
Monitoring Camera (OMC), to provide photometry in the 
Johnson V-band simultaneously to the high-energy observations 
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( |Mas-Hesse et al.||2003| l. In addition to monitoring the primary 
high-energy targets, OMC has been observing serendipitously 
tens of thousands of potentially variable objects within its field 
of view, preselected from existing catalogues of variable stars. 

In this paper we present the first INTEGRAL-OMC cata- 
logue of optically variable sources (OMC-VAR hereafter) that 
contains photometric data in the V-band, as well as variability 
and periodicity data for more than 5000 objects observed from 
October 2002 to February 2010. The catalogue will be updated 
as new data are being analysed. At the end of the mission we ex- 
pect to produce the final catalogue with variability information 
for more than 25 000 objects monitored over a long time period, 
with consistent and well calibrated photometry from space. 



2. The Optical Monitoring Camera onboard 
INTEGRAL 

The Optical Monitoring Camera (OMC) is a small refractive 
telescope with an aperture of 50 mm focused onto a large for- 
mat CCD (1024 X 2048 pixels) working in frame transfer mode 
(1024 X 1024 pixels imaging area). With a field of view of 5°x5° 
OMCprovides photometry in the Johnson V-band (centred at 
5500 A) and for sources from V ~ 7 mag (for brighter sources 
saturation effects appear) to V ~ 16-17 mag (limit magnitude 
for 3cr source detection and 200 s integration time). Typical ob- 
servations are performed as a sequence of different integration 
times, allowing for photometric uncertainties below 0. 1 mag for 
objects with V < 16. The accuracy is limited to 0.01 mag by 
the flat field correction matrix. OMC was designed to observe 
the optical emission from the prime targets of the gamma-ray in- 
struments onboard INTEGRAL and has the same field of view 
(FoV) as the fuUy coded Fo V of the INTEGRAL X-ray monitor 
(JEM-X; [Lund et al.|2003| , and it is co-aligned with the central 
part of the gamma-ray instruments larger fields of view, the IBIS 
imager (Ubertini et al. 2003) and the SPI spectrometer (jVedrenne 
et al.|2003[ ). The need to have a FoV as large as possible led to 
a rather large angular pixel size (17'.'5 x 17'.'5), with a Gaussian 
Point Spread Function with FWHM ~ 1 .4 pix. 



2.1. The OMC Input Catalogue 

In addition to monitoring the INTEGRAL primary high-energy 
targets, OMC has the capability to observe serendipitously 
around 100 sources within its FoV. These targets have to be pre- 
selected on ground, since the limited telemetry available to OMC 
(only ~ 2.2 kbps) does not allow to download full images to 
ground. These additional targets are automatically selected from 
the OMC Input Catalogue (OMC-IC hereafter) by a specific 
OMC Pointing Software (OMCPS) running at the INTEGRAL 
Science Operation Centre (ISOC). Only 11x11 pix windows of 
the CCD containing those objects are transmitted to ground. The 
OMC-IC was compiled before the start of the mission, and is 
described in detail in D omingo et al.| ( i2003) . It has been updated 
several times to include new targets of interest or to improve the 
accuracy in coordinates, especially of the high-energy objects. 

The current version of the OMC-IC contains over 500000 
entries, namely: 

- astrometric and photometric reference stars, 

- known optical counterparts of gamma-ray sources, 

- known optical counterparts of X-ray sources, 

- point-like X-ray sources detected and catalogued by ROSAT, 

- quasars observable by the OMC, 



- known additional AGNs and 

- a big fraction of the known variable stars (including eruptive 
variable stars, novae and cataclysmic stars). 

The OMC-IC includes both known variable stars as well as 
suspected variable stars. Special care was devoted to the com- 
pilation of extragalactic objects (AGN, QSO, starburst galaxies, 
narrow emission line galaxies, etc.) as they are potential high- 
energy emitters, and in many cases show significant variations 
of their optical flux. Indeed, Beckmann et al. ( |2009) compiled 
the second INTEGRAL AGN catalogue including optical infor- 
mation from OMC for 57 AGN. 

In the compilation of optical sources, mostly available before 
the INTEGRAL launch, the SIMBAD database and the follow- 
ing catalogues were used: 

1. Combined General Catalogue of Variable Stars, 4.1 (GCVS; 
Kholopovetal.|I998[ ). 



It contains a catalogue of all known Galactic variable stars 
prior to 1997, a catalogue of extragalactic variable stars, a 
catalogue of supernovae and a catalogue of suspected vari- 
able stars not designated as variables prior to 1980 (pub- 
lished in The New Catalogue of Suspected Variable Stars, 
Kukarkin & Kholopov|I982] l. 



New Catalogue of Suspected Variable Stars. Supplement 
( |Kazarovets et al.|I998l ). 

It contains 1 1 206 stars suspected of variability which were 
not designated as variables prior to 1997. 
The 74th Special Name-list of Variable Stars ( |Kazarovets| 
et al. 1999). 

It contains 3157 variable stars whose variability was dis- 
covered by the Hipparcos mission. All the stars satisfy the 
GCVS variability criteria. 

The Hipparcos Catalogue ( Ferryman et al.|1997 1. 

The 1 1 597 variable stars in the Hipparcos catalogue were 

included. 

Variable stars in the Tycho photometric observations. 
To search for variability among faint stars of the Tycho cata- 
logue ( |Hoeg et al.|I997| l, |Piquard et al.| ( |200 1) made a treat- 
ment that took into account truncated detections and cen- 
sored measurements. The list contains 1091 stars suspected 
to be variable stars. 



7. 



Quasars and Active Galactic Nuclei, 8th Ed. ( Veron-Cetty & 
|Veron|I998| l. 

The Active Galactic Nuclei Catalogue (Padovani priv. 
comm.). 

This catalogue includes 12021 quasars and active galaxies 
and is heavily based on the catalogue of Quasars and Active 
Galactic Nuclei, 7th Ed. ( [Veron-Cetty & Veron|19 96). It also 
includes the BL Lac Catalogue ( [Padovani & Giommi^l995] l 
updated with BL Lac's discovered in 1996, and the radio 
galaxies in the 1 Jy, S4, and S5 radio catalogues. The AGN 
catalogue reports V magnitudes almost for all of their ob- 
jects. Nevertheless, in some cases V magnitudes were de- 
rived from B magnitudes by assuming a {B - V) colour in- 
dex. 

Narrow Emission Line Galaxies (Kunth priv. comm.). 
Compilation of 441 objects dominated by intense starburst 
activity. These objects show strong optical emission lines 
and/or very blue stellar continuum, both being tracers of 
young, massive stellar populations and therefore candidates 
for host supernovae. 

Catalogue of Cataclysmic Variables ( jDownes et al.||2()0I 



This catalogue contains 1134 sources that were included in 
the OMC-IC. 
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Fig. 1. Visual description of the process of period determination for two different cases: a standard situation (analysis of lOMC 
0460000022) on the left and an example of aliasing (lOMC 5228000021) on the right. In both cases, the plots ai-e: Top: The light 
curve without folding. For the light curve affected by aliasing the large temporal gap can be appreciated. Middle: Periodogram from 
to 10 days. The peaks are much narrower in the standard case. In the aliasing example, the smoothed periodogram taking into 
account the aliasing effect (see text) is overplotted in green. Bottom: Detail of the method used to determine the uncertainty of the 
period. In red, the continuum plus gaussian fit to the peak is plotted. 



10. Candidates for classical T-Tauri stars 

Classical T Tauri stars are good targets for OMC because 
they emit in the optical band in a periodic or non periodic 
way, but always have some type of variability related to their 
physical conditions. The amplitudes of the light curves are 
in the range 0.01-3 mag. Periods are often ~1 day, so the 
monitorization from a space platform should overcome the 
lack of observing them from the ground. 732 candidates for 
classical T-Tauri stars were selected ( |Caballero-Garcia|2003| l 
based on the segregation shown by these stars in a 2MASS 
(/ - H) vs. {H - Ks) colour-colour diagram ( |Lee & Chen| 
[2002) 1. 

1 1 . Candidates for cool dwarfs ofG-K-M spectral type. 

It is broadly known that cool dwarfs of G-K-M spectral type 
can be variable at a certain degree (mainly in the range 0.05- 
0.1 mag. and not greater than 0.2 mag.), showing very inter- 
esting properties of the behaviour of this kind of stars. A se- 



lection of 35 101 candidates was made by Caballero-Garcia 



( 2003 | l based on proper motions and the 2MASS (7 - 
{H - Ks) colours. 



H) and 



3. Data analysis 

After the proprietary period of one year, all INTEGRAL data 
are open to the scientific community. At the moment of writing, 
INTEGRAL has been in orbit for more than nine years and the 
OMC database {http://sdc.cab.inta-csic.es/omc/\ Gutierrez et al. 
20041 now contains observations for more than 130000 objects, 
of which more than 70 000 have light curves with at least 50 pho- 
tometric data points. The data included in this catalogue contain 
observations from the beginning of the mission until February 
2010. 



3. 1 . Selection of the sources and cleaning of the light curves 

For this work, sources with more than 300 photometric points 
have been selected from the OMC database, in order to deal 
with light curves with enough points to study their variability 
and when possible, their periodicity. With the aim of improving 
the quality of the data some selection criteria have been applied 
to the individual photometric points. 

First of all, we have checked if the coordinates in the OMC- 
IC agree with the coordinates given in SIMBAD. The OMC 
database allows the user to select between two different centroid 
methods to extract the photometry, which corresponds to dif- 
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Fig. 2. Histogram of magnitudes, covering the range from 7.10 
to 16.27, with a peak of the distribution at around 12.3 mag. In 
this and subsequent histograms the fraction of sources for which 
we have determined a period in this work are represented with 
dark green. 



ferent options in the standard OMC Off-line Science Analysis 
(OSA) softwar e distributed by the INT EGRAL Science Data 
Centre (ISDC; [Courvoisier et al. | |2003p . If Source coordinates 
(default) is selected as the centroid method, the pipeline looks 
for the centroid allowing a maximum offset of 10 arcsec with 
respect to the given coordinates. If the centroiding algorithm 
fails by any reason, the extraction mask is nevertheless re- 
centered right on the source coordinates as given in the OMC- 
IC. However, if brightest pixel is selected as the centroid method, 
the position of the extraction mask is re-centered on the brightest 
pixel within the central 5x5 pixel section. For the sources with a 
difference between the coordinates of the input catalogue and the 
coordinates given by SIMBAD greater than 10 arcseconds, the 
brightest pixel method has been used. Moreover, in all the cases, 
those photometric points with a distance to the coordinates of 
SIMBAD greater than 10 arcseconds have been rejected. 

Secondly, we have looked for duplicated entries in the cat- 
alogue. There are sources in the input catalogue with different 
OMC identifiers, but with the same SIMBAD name or very close 
coordinates. These two situations have been taken into account 
and the duplicated sources have been removed from the list. 

Typical OMC integration times are 10 s, 50 s and 200 s, al- 
though at the begining of the mission shots of 10 s, 30 s and 
100 s were obtained. Therefore, photometric points with differ- 
ent integration times can be found in the database. The longest 
integration times led to saturation of the brightest objects. For 
this reason, for each light curve we have rejected the photomet- 
ric points that could be saturated. In general, all the photometric 
points within light curves where the median brightness is less 
than 7.1 were rejected. When the median brightness was less 
than 8.3, those photometric points with exposure times greater 
than 20 s were eliminated. If the median brightness was between 
8.3 and 8.8, only the photometric points with exposure times of 



10 or 30 s were selected. For light curves with median bright- 
ness between 8.8 and 9.6 the exposures of 100 and 200 s were 
rejected. When the median brightness was between 9.6 and 10.3, 
the photometric points with exposure times of 200 s were re- 
moved. And finally, for those sources with median brightness 
greater than 10.3 we considered that no problems of saturation 
should appear and all photometric points were considered. 

In a similar way, some detection limits were defined to avoid 
photometric points with too low signal-to-noise ratio (SNR). For 
each light curve, those photometric points with exposure times 
of 10 s and a median brightness greater than 12, exposure times 
of 30 s and a median brightness greater than 13 and exposure 
times of 50 s and a median brightness greater than 14 have been 
removed. 

Furthermore, measurements with low SNR were excluded. A 
minimum SNR of 10 was required for the shortest integrations 
(10, 30 and 50 s) while for the longest ones (100 s and 200 s) a 
criteria of SNR > 3 was imposed. 

The observations were filtered from the effects of cosmic 
rays hits and some occasional detector artifacts. Each photomet- 
ric point was compared with the 20 closest ones within 3 days 
(to avoid the timing gaps). For these 20 points, the standard devi- 
ation was computed and a 5-sigma rejection criteria was applied. 
We have checked manually that this algorithm removed most of 
the outliers. Finally, the photometric points that contained any 
hot pixel within their extracting mask of 3 x 3 pixels were re- 
moved. 

In this first version of the OMC catalogue we have anal- 
ysed only those objects with a given SIMBAD object type in 
the OMC-IC. Applying these filters we obtained 6071 sources 
with high-quality light curves to be studied. In the next version 
of the OMC catalogue we will include all objects in the OMC 
database complying with the previous requirements, indepen- 
dently of their assigned type, estimated to be around 25 000. 



3.2. Detection of variability 

To detect which light curves showed variability we have used 
the goodness-of-fit test. This test compares a null hypothe- 
sis to an alternative hypothesis. If the value obtained for the test 
statistic is greater than a value corresponding to the chosen sig- 
nificance level of the distribution with k = n-1 degrees of 
freedom, being n the number of data points, then we reject the 
null hypothesis. We have used for each source the null hypothe- 
sis of being non-variable and having a constant magnitude equal 
to the mean. The alternative hypothesis is that the star is variable. 
Then, we calculated the x^ and the significance a. This signif- 
icance gives the probability of being wrong when rejecting the 
null hypothesis (the source is constant). We have considered as 
variable those sources with a < 0.05 (a lower limit probability of 
being variable of 95%). Following this criterion we have identi- 
fied 5263 sources showing variability in their OMC light curves, 
which constitute our present catalogue. Though all 6071 sources 
were expected to be variable, we have not identified variability 
for 808 of them, according to the criteria we have applied. In 
some cases this can be due to the time span of the OMC ob- 
servations, and in other cases the amplitude is likely below the 
OMC variability detection threshold used for this first version of 
the catalogue. 
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Fig. 3. Histogram of amplitudes (calculated as the difference be- 
tween the magnitude at minimum and the magnitude at maxi- 
mum of the light curve) with a peak of the distribution at around 
0.84 mag. In dark green are represented the sources for which 
we have determined a period in this work. 

3.3. Study of the periodicity 

To determine which sources are periodic and to derive their pe- 
riods, an algorithm based on the Phase Dispersion Minimization 

1978| l has been developed. 



(PDM) technique (StelUngwerf 



Given the variety of variability patterns present in our data we 
decided not to use the Lomb-Scargle (L-S) periodogram tech- 
nique by 'Scarglel ( |1982| ). This method evaluates the discrete 
Fourier transform for nonevenly sampled data. Since this is 
equivalent to assuming a prior sinusoidal light curve, we consid- 
ered it not adequate for most eclipsing and non-sinusoidal light 
curves. 

We explored a wide range of frequencies corresponding to 
periods from 0.8 h (a lower limit of the periods that can be de- 
tected by OMC, considering the typical exposure times), up to 
half of the time range of the observations (the maximum period 
to have at least two cycles in the folded light curve). We first used 
a frequency increment of j^, being B the total time baseline of 
the observations, and in a second iteration around the first trial 
period, we used a period increment of 2P^/{nB) (see Gimenez 
& Garcfal|1983 1, where P is the trial period and n is the num- 
ber of points in the light curve. For each period, the time-folded 
light curve is divided into a series of bins and the variance of 
the amplitude within each bin with respect to a mean curve is 
computed. In our algorithm the number of bins is variable and 
it depends on the quality of each light curve. For each bin, the 
mean value is calculated and linear interpolations are done be- 
tween the means in order to obtain one average curve for each 
folded light curve. The ratio between the sum of the bin vari- 
ances and the overall variance of the data set is called 0, and the 
period that minimizes this parameter has been considered to be 
the best estimate. 

With the values of the trial periods and their corresponding 
values of 0, we computed the periodograms. The uncertainty in 
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Fig. 4. Histogram of the periods derived. Typical values vary 
between a few hours and 10 days, with a peak of frequency at 
around 15 h. 
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Fig. 5. Comparison of periods derived from OMC data (this 
work) with the values compiled by the VSX. The solid line rep- 
resents PoMC - Pvsx while the dashed line marks the condition 
PoMC - ^Pvsx- 



the period was calculated by fitting a gaussian to the peak cor- 
responding to the minimum value of in the periodogram. The 
1-sigma error in the determination of the centroid of the gaus- 
sian is provided for each source as the standard uncertainty of 
the period and corresponds to a confidence interval of 68.3% 
(see Fig. [TJ. For most of the sources the process above works 
fine, but in some cases aliasing becomes a problem. Aliasing 
occurs when a curve is observed with an insufficient sampling. 
For some sources in the catalogue there are temporal gaps in the 
observations much greater than the value of the period, and the 
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Table 1. Contents of the OMC-VAR catalogue, accessible at http://sdc.cab.inta-csic.es/omc/var/ and at CDS services. 



Field 


Source 


Description 


lOMC 


OMC 


OMC Identifier 


SIMBAD name 


SIMBAD 


SIMBAD Basic Identifier 


RA 


SIMBAD 


Right Ascension in sexagesimal units (J2000) 


DE 


SIMBAD 


Declination in sexagesimal units (J2000) 


RAdeg 


SIMBAD 


Right Ascension in degrees (J2000) 


DEdeg 


SIMBAD 


Declination (J2000) in degrees (J2000) 


otype 


SIMBAD 


SIMBAD object type 


sptype 


SIMBAD 


SIMBAD spectral type 


vartype 


VSX 


VSX variability type 


vargroup 


VSX, SIMBAD, others 


General Variability Group 


V mag 


OMC 


Median V magnitude 


magerr 


OMC 


Mean of the magnitude error 


magMax 


OMC 


Magnitude at maximum brightness 


magMin 


OMC 


Magnitude at minimum brightness 


per 


OMC 


Period of variability in days 


errper 


OMC 


Uncertainty of the period in days 


timebase 


OMC 


Time base of the observations in days 


tagcont 


UCAC3, NOMADl, SIMBAD, VSX 


Tag of contamination 



time intervals of the different observations are not long enough 
to determine the period accurately. For a given estimated period 
Pi, if the temporal gap, fgap (the time when the source has not 
been observed), is such as fc,ap - N x Pi (with big enough), 
there will be many significant peaks in the periodogram, close to 
the deepest one, corresponding to that values of the period that 
make fg^p = (A^ + 1) x or fgap = (A^ - 1) x (see Fig.[T] 
right panel). In these cases we have smoothed the periodogram 
with a window size wider than the distance between these false 
peaks. This smoothed curve is used to determine the limit of sig- 
nificance of the peak and the half width of this peak is used as 
initial value of sigma in the gaussian fit. Finally, a visual inspec- 
tion of each folded light curve was performed to identify and 
correct potential problems due to ambiguities in the results, such 
as deriving the half or double period in some eclipsing binaries 
or pulsating stars with symmetrical light curves. 

The method used in this work yields only the primary peri- 
ods. To obtain second or multiple periods, if present, a manual 
analysis of each source is required and this is out of the scope of 
this paper 



veloped three different methods. The first one was to take into 
account the photometric information of the sources close to the 
extraction aperture present in the UCAC3 (Zacharias et al. 2010| l 
and NOMADl (Zacharias et al. 20041 catalogues. The second 



We show in Fig. A. 5 the light curve of Her X-1, one of 
the brightest and most studied persistent X-ray binary pulsars. 
The system displays a great variety of phenomena at different 
timescales, including pulsations at 1.24 s, eclipses at the orbital 
period of 1.7 days, and a 35 day cycle in the X-ray intensity 
( Leahy & Dupuis|2010| l. Despite the large gap between the two 
epochs at which this object has been observed by OMC, our 
method retrieves properly the orbital period. Nevertheless, the 
length of the gap, compared to the short orbital period, implies 
that the period has been derived with a rather low accuracy. 

3.4. Photometric contamination by nearby stars 

For the data analysed in this study, fluxes and magnitudes were 
derived from a photometric aperture of 3x3 pixels (1 pixel - 
17'.'504), slightly circularized, i.e., removing | of a pixel from 
each comer (standard output from OSA). Therefore the com- 
puted values include the contributions by any other source in- 
side the photometric aperture. We flagged in a column of the 
catalogue those sources that might be affected by this kind of 
photometric contamination. Due to the variable nature of the 
sources in the catalogue and the difficulty to accurately estimate 
the contamination by the nearby and background stars, we de- 



one was to compare the V magnitude provided by SIMBAD (for 
variable sources it is usually the magnitude at maximum bright- 
ness) with the lowest magnitude we measure for each source, 
and the third one was to compare the median magnitude in VSX 
with the median magnitude we measure. We considered a source 
significantly contaminated when the measured flux was a fac- 
tor of 2 higher than expected. This is equivalent to a difference 
between expected and observed magnitude of more than 0.75 
mag. While this criterion is not very restrictive in terms of pho- 
tometric accuracy, it allows to study variability and search for 
periodicity for objects with a significant degree of contamina- 
tion. We want to stress nevertheless that this procedure has the 
caveat that in some cases the variability/periodicity found might 
be originated by the contaminating object(s), and not by the tar- 
get source. One of the objects mostly affected by blending is HH 
Nor, shown in Fig. A. 6 whose variability pattern is completely 
dominated by the variability of the brighter, nearby object VSX 
J 154329.4-5 15037, which is indeed a variable RR Lyr star with 
a much shorter period. Being located at just 12", a fraction of 
an OMC pixel, the intrinsic variability of this contaminating star 
completely hides the variability of the target source, HH Nor 

The results of these three methods are coded in the column of 
contamination with three numbers (a reference to each method) 
where 1 means that the source is not contaminated (under the 
conditions defined above), 2 means that it is contaminated and 
3 means no photometric information has been found for that 
source. Because of the inherent difficulties of estimating the con- 
tamination of variable sources, this flag is for guidance only and 
to study individual sources a more exhaustive analysis of the sky 
field around each source is recommended. 

4. Contents of the OMC-VAR catalogue 

Table [T| describes the contents of the OMC-VAR catalogue, 
which is fully accessible online. In this first release we provide 
information about the variability of 5263 objects distributed all 
over the sky (see Fig.|6]l. The majority of the objects are located 
in the Galactic plane and Galactic bulge, as well as on specific 
areas of the sky. This distribution reflects the integration time 
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Fig. 6. Distribution in galactic coordinates of all the sources in the catalogue. The green crosses represent pulsating stars, the red filled points 
correspond to eclipsing binaries, the brown plusses are eruptive stars, the orange complex plusses represent rotating stars, the inverted dark blue 
filled triangles represent cataclysmic variables, the light blue filled stars are X-ray binaries, the yellow filled points correspond to objects simply 
classified as variable stars, the purple filled pentagons represent extragalactic objects and the empty purple squares are other types of objects. 



density maps of INTEGRAL, driven by the location of the most 
interesting high-energy sources. 

For each source we include the median value of the V mag- 
nitude, the mean of the photometric errors of the points and the 
magnitudes at maximum and minimum brightness, estimated as 
the 2nd and 98th percentiles, respectively, of the points in the 
light curve. For those sources classified as periodic or probably 
periodic, we provide the best estimate of the period and its cal- 
culated uncertainty. 

The values of the median magnitudes go from 7.10 mag to 
16.27 mag. In Fig. [2] we show the histogram of the median mag- 
nitudes, showing a peak at around 12.3 mag. A similar plot for 
the amplitudes (calculated as the difference between the mag- 
nitude at minimum and the magnitude at maximum of the light 
curve) is shown in Fig. [3] The amplitudes go from 0.06 mag to 
7.41 mag with a mean value of 0.84 mag. 

The time base for each light curve is also provided. These 
time bases have a minimum value of 0.69 d, a maximum value 
of 2815.19 d and a mean value of 1862.0 d. 

We have determined periods for 1337 sources, out of the 
5263 variable objects compiled in the catalogue. The distribu- 
tion of periods is shown in Fig.|4] Typical values vary between a 
few hours and 10 days, with a peak of frequency at around 15 h. 
Some sources present periods of up to some hundreds of days, 
but there are not many of them given the time span and distri- 
bution of the observations. The VSX database also provides in- 
formation about the periodicity of the sources and we have used 
this information to compare with our results. We have derived 
periods for 175 objects whose periodicity was unknown and in 
many other cases, we have improved the results with respect to 
those found in the VSX database, as can be seen in Appendix A. 

In Fig. |5] we compare the periods derived from OMC data 
in this work with the values compiled in the VSX, when both 
were available. For the large majority of objects the values are 
consistent, though we want to stress that in many cases, even a 
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Fig. 7. Histogram of the different groups of variability types 
present in the catalogue. There are a big number of objects which 
are unclassified, or just classified as variable objects. 

very small improvement on the determination of the period is 
enough to yield a satisfactory folding of the light curve (see the 
examples in Appendix A). A dashed line marks the points for 
which PoMC = 2 X Pysx^ to identify those cases of binaries with 
eclipses of similar shape that could have been misidentified in 
the VSX. In the rest of the cases we have verified that the period 
given in the VSX was not consistent with the OMC light curves. 
This happened sometimes for systems with multiple periods due 
to different variability mechanisms like rotation, orbital. X-rays, 
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Fig. 8. Histogram of the subgroups of variability types present 
in the catalogue. Top: Histogram of the different kinds of pulsat- 
ing stars. Bottom: Histogram of the subtypes of eclipsing binary 
stars. 



etc. One example is Her X-1 for which, as we mentioned in 
section r 



3.3 the period given by VSX does not correspond to the 



orbital one to which OMC is more sensitive (see Fig. A. 5 1 



4.1. Classification 

A big variety of variable objects have been observed by OMC 
and are present in this version of the catalogue. In Fig.|6]the dis- 
tribution in galactic coordinates of the different kind of objects 
is plotted. In order to do a statistical study of the type of objects 
in the catalogue we have used the information available in the 
literature. When available, the object classification has been ex- 
tracted from the VSX. Otherwise, the type of variability has been 
taken from the SIMBAD database. According to the above clas- 
sifications we have defined different variability groups and have 
assigned each source to one of them. In the case of pulsating 
stars we have defined 10 additional subgroups, and in the case 
of eclipsing binaries we have included 4 additional subgroups, 
based only on the shape of the light curve: 

- 2631 pulsating stars (Pu), including 



• 297 Cepheids (Ce), 

• 261 RR Lyrae (RR) pulsators, 

• 31 6 Scuti pulsators (6 Set), 

• 32 RV Tau stars (RV), 

• 28/3 Cephei stars (j3 Cep), 

• 26 SPB pulsators (SPB), 

• 680 Mira type stars (M), 

• 591 semi -regular stars (SR), 

• 630 slow pulsators (L), 

• 49 simply classified as pulsators (other). 

- 11 32 eclipsing binaries (EB), including 

• 728 Algol (J3 Persei)-type eclipsing systems (EA), 

• 191 /3 Lyrae-type eclipsing systems (EB), 

• 1 17 W Ursae Majoris-type eclipsing variables (EW), 

• 96 simply classified as eclipsing (E) 

- 53 X-ray binaries (XB), 

- 78 cataclysmic variables (CV), 

- 312 eruptive stars (Er), 

- 144 rotating stars (Ro), 

- 79 Be stars (Be), 

- 579 stars just classified as variable (V*), 

- 80 extragalactic sources (G), and 

- 162 objects of other types. 

We show in Fig. |7]the histogram corresponding to the dif- 
ferent groups, clearly dominated by the groups of pulsating stars 
and eclipsing binaries. In Fig.[8]we detail the histograms for the 
subgroups of pulsating stars and eclipsing binary systems, the 
latter dominated by Algol-type systems. 

Some examples of the kind of objects that can be found in 
the catalogue are shown in Appendix A, where we have plot- 
ted some representative charts including the DSS image around 
the target, and their light curves, unfolded and folded with the 
periods we have derived and/or with the catalogued ones. 

5. Summary 

We have analysed the optical photometric light curves, covering 
the time from the beginning of the mission to February 2010, of 
the 6071 objects present in the OMC database that satisfy the 
criteria explained in section [3] with the aim of characterizing 
their expected variability and potential periodicity. 

In this first version of the INTEGRAL-OMC catalogue of 
optically variable sources we have detected variability for 5263 
objects, with median brightnesss in the range from 7.1 to 16.27 
in V. 

For 1337 of them we have been able to derive a period with 
enough significance. The periods found vary from a few hours to 
some hundreds of days, with a peak in the histogram of typical 
periods around 15 h. We inspected visually the light curves to 
check the validity of the estimated periods, rejecting those false 
positives yielded by the PDM algorithm. 

For each source in the catalogue we include the median value 
of the V magnitude, the mean of the photometric errors of the 
points and the minimum and maximum value of V, estimated as 
the 98th and 2nd percentiles, respectively, of the points in the 
light curve. For those sources classified as periodic or probably 
periodic, we provide the best estimate of the period and its calcu- 
lated uncertainty. We have also included a set of tags to identify 
those sources whose photometry can be significantly affected by 
nearby stars, as well as some ancillary data. 

This catalogue contains data useful for the study of several 
classes of variable stars: eclipsing binaries, pulsating stars, ac- 
tive stars, cataclysmic variables. X-ray binaries, etc. The most 
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frequent objects in the present compilation are pulsating stars 
and eclipsing binary systems, according to the classification pro- 
vided by the VSX catalogue and SIMBAD. A deeper study of 
some of the sources presented here is being carried out and wiU 
be presented in forthcoming papers. 

The next version of the INTEGRAL-OMC catalogue of op- 
tically variable sources will include all other objects in the OMC 
database complying with the same photometric requirements, es- 
timated to be around 25 000. 
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Appendix A: Charts provided for each object in the 
catalogue 

Some illustrative examples of the charts provided with the cata- 
logue for different kinds of objects. The information about type 
of object, variabihty type and spectral type on the top of each 
plot comes from Simbad and VSX. In all cases we show the DSS 
red image of the field of view in the top left panel (red circles 
represent the OMC photometric aperture), and the OMC V-band 
light curve without folding in the top right one. The time is the 
Barycentric INTEGRAL Julian Date expressed in Barycentric 
Dynamical Time (TDB). To convert this time to Barycentric 
Julian Date expressed in TDB, you have to add 2,451,544.5 d. 
Similar charts and the complete light ciu'ves in machine readable 
format for all the sources contained in the OMC-VAR Catalogue 
can be retrieved from 
http://sdc. cab. inta-csic. es/omc/. 
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Fig. A.l. IT Nor, an Algol type eclipsing binary. When not indicated otherwise, we show in the different charts the OMC light curve 
folded with the period estimated in this work (bottom left), and the OMC light curve folded with the period given by the VSX 
catalogue (bottom right panel). In this case the period derived from the OMC light curve improves significantly the value provided 
by the VSX and yields a satisfactory folding. 
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Fig. A.2. For some objects, Uke the RR Lyr star S Com, a very small increase in the accuracy of the period is enough to get a perfect 
folding of the light curve. 
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Fig. A.3. For a large fraction of objects, the period derived from OMC data and that compiled by VSX yield similar results when 
folding the light curve. 
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Fig. A.4. New period determination from this work for V65 1 Ori. Moreover, the variability classification in VSX should be corrected 
from RR Lyrae to eclipsing binary, according to our light curve. 
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Fig. A.5. In the case of Her X-1, the VSX period corresponds to the 35 day cycle in the X-ray intensity, while the OMC optical 
light curve is dominated by the orbital eclipses. 
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Fig.A.6. The photometry of the intended target, HH Nor, is heavily dominated by the brighter, nearby object VSX J 154329. 4— 
515037, which is indeed a variable RR Lyr star (RRAB) with a period of 0.59 days. In fact, HH Nor is an eclipsing binary with 
mean V magnitude of 10.9 and a period of 8.58 days. 
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Fig.A.7. In the light curve of this high-mass X-Ray binary we see that a long term variation during the nearly 3000 days of 
monitoring dominates over the orbital period given by VSX. When there is no determination of the period from this work, the chart 
of the corresponding object includes a histogram of the observed magnitudes. The shape of these histograms helps to understand 
the variability pattern of the different objects. 



lOMC 4368000068 8C 0716+714 

Type: BLL Var Type: BLLAC SP Type: 

= 14.3 mag <err_V> = 0.1 mag 



12,5 

13.0 

"Si 13.5 
o 

> 14.0 
14.5 
15.0 



1400 



1450 1500 1550 
Borytime (do/s) 



1600 



300 
250 
2DD 
150 
100 
50 




13.0 13.5 



14.0 14.5 
V (maq) 



15.0 



Fig. A.8. The OMC light curve of this blazar shows two different states, separated by around 150 days. They can be clearly observed 
in the magnitude histogram. 
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